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We give an algorithm allowing to construct bases of local unitary invariants of pure 
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O 
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—^ ' fc-qubit states from the knowledge of polynomial covariants of the group of invertible 

f~^ , local filtering operations. The simplest invariants obtained in this way are explicited and 

^^3 ' compared to various known entanglement measures. Complete sets of generators are 

^p , obtained for up to four qubits, and the structure of the invariant algebras is discussed in 

detail. 

-(— > 
^ . 

qh, 1. Introduction 

From a mathematical point of view, Quantum Information Theory deals with finite di- 
^^ . mensional Hilbert spaces, the state spaces of finite fc-partite systems, which have the 

H ' special form 

H = Vi®V2®---®Vk, (1) 

where Vi is the finite dimensional state space of the ith part (or particle) of the system, 
most of the time assumed to be a qubit, which means that dim Vi — 2. 

The interesting non-classical behaviors on which the theory is based already occur for 
two-qubit systems, for the so-called entangled states, those ij} ^Vi®V2 which cannot be 
written in the form wi (g)W2 • The properties of such states are the basis of the EPR paradox 
IjEinstein et al. 1985|l . and since its discovery, the entanglement phenomenon has been 
thoroughly investigated by physicists, cf. IjBell 1966IIClauser et al. 1969l|Aspect et al. 1982| 

IBennett and Wiesner 1992)l . and more recently by mathematicians, e.g., ( |Brylinski and Brylinsky 2002| 

[El^chko 20021 peyer and Wallach 2002| ). 

There is, however, no general agreement on the definition of entanglement for systems 
with more than two parts. Klyachko has proposed ( |Klyachko 2002|rKlyachko and Shumovsky 2003| ) 
to regard as entangled the states which are semi-stable for the action of the group of in- 
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vertible local filtering operations, also called SLOCC^, 

G = SL{Vi) X • • • X SL{Vk) (2) 

in the sense of geometric invariant theory, which means those states on which at least one 
non trivial G-invariant polynomial does not vanish. The point of introducing geometric in- 
variant theory is that this theory provides methods for characterizing such states without 
explicitly computing the invariants. In order to explore the significance of this property, 
the invariants have been explicited in the simplest cases (up to 4 qubits, and 3 qutrits, 
with partial results for 5 qubits, (Luque and Thibon 2003^|Luque and Thibon 2005|IVerstraete et al. 2002|l ). 

One would also like to quantify entanglement. The non- locality properties of an entan- 
gled state does not change under unitary operations acting independently on each of its 
sub-systems. The idea of describing entanglement by means of local unitary invariants is 
explored in l|Grassl et al. 19981 . see also (Schlienz and Mahler 1996: Schlienz and Mahler 1995|l . 
However, except for the simplest systems, there are far too many orbits and a complete 
classification is out of reach. 

An intermediate possibility is to look at the G-orbits. The knowledge of the G-invariant 
polynomials is not sufficient to separate the G-orbits, and in general, one has to look for 
the covariants in the sense of classical invariant theory. 

In HBriand et al. 2003|l . the algebra of G-covariants for 4 qubits has been investigated, 
and a complete set of generators has been obtained. 

In the present article, we will explain how these results can be applied to the calculation 
of bases of unitary invariants. As an application, we compute bases of the spaces of local 
unitary and special unitary invariants of degree 4 of A: qubits for arbitrary fc, and recover 
the results of l|Grassl et al. 2002)l for 3 and 4 qubits. 

The paper is organized as follows. In Section Owe recall some background on SLOCC 
covariants, and we describe a method allowing to obtain from them local unitary (LUT) 
and special unitary (LSUT) invariants. Sectionals devoted to the computation of the 
simplest LUT and LSUT invariants from SLOCC covariants. Finally, we give some ex- 
amples and applications in Section [S] 



2. Invariants and covariants of qubit systems 

2.1. Group actions on state spaces 

Let y = C ^ be the local Hilbert space of a two-state system (a qubit), and T-i = V®^ be 
the state space of a system of k qubits. We shall regard it as the natural representation of 
the group G = Gslocc = SL{2, C )'^, known in quantum information theory as the group 
of reversible local filtering operations, or stochastic local quantum operations assisted by 
classical communication (SLOCC) IjBennett and Wiesner 19921 IDtir et al. 200 l|l . This is 
a semisimple complex Lie group, whose representation theory follows immediately from 
that of SL{2,C). The maximal compact subgroup oi G is K = Glsut = SU{2)^, the 

t for Stochastic Local Operations assisted with Classical Communication. 
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group of local special unitary transformations. We shall also be interested in arbitrary 
local unitary transformations, which form the group U = Glut — U{2)''. 
If |j)i J = 0, 1 is a basis of V, a state |^) can be written as 

1 
I*) = X! ail^2■■■ik\^l^2■■■ik) (3) 

il,...,ik=0 

where, as customary in the physics literature, 

\iii2---ik) = \h) ® •••«' \ik) ■ (4) 

It will be convenient to interpret such a state as a multilinear form 

/(x) = /(2;W,...,x('=))= J2 «^.^.-^.4'^---^i^ (5) 



n,...,ifc=0 



where a;^-'^ = (xq ,Xi ) are pairs of variables. 

The action of a fc-tuple of matrices g = (g*-^\ . . . ,5*-'^'') on the various vector spaces 
introduced so far is defined by gx = x', x'^'^ = g^^^'x^^^ and the components (1^^^^...^^^ of 
/' = g/ are defined by the condition 

ii,...,ifc ii,...,ifc 

In the following, we shall be interested in LUT and LSUT invariants of a state |'0), i.e. 
polynomial functions /(a, a) in the components of jf/'), such that 

/(a,a)-/(a',a') (7) 

where a[ ...^ are the components of /' — g/ for g a LUT or a LSUT. Our main point 
will be the application of the SLOCC invariant theory to the calculation of such unitary 
invariants. 



2.2. SLOCC Invariants 

The SLOCC invariants are the holomorphic polynomials /(a) such that /(a) = /(a') 
for g G GsLocc- Of course, the squared modulus |/p of a SLOCC invariant is a LSUT 
invariant, but only a small subset of unitary invariant are of this form. 

The methods of classical invariant theory can be applied to the determination of the 
SLOCC invariants of k qubits for small k. An important preliminary step is the deter- 
mination of the Hubert series 

h{t) =^t'^ dim 5'^(?^)^«^occ ^ (g) 

the generating series of the dimension of the space of homogeneous polynomial invariants 
of degree d. 

For A; = 3, the only fundamental polynomial invariant of three qubits is known since 
the nineteenth century, it is the Cay ley hyperdeterminant ( |Le Paige 1881| ), see also 
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( IMiyake 2003| ): 

■pv J./4-1 (2 2,2 2,2 2,2 2\ 

ueti^^j — (.agpoam + aQoia^^Q + aoioflioi + Qoii^iooy 

— 2(aoooaooiaiioaiii + aoooaoioQ-ioiaiii + aoooaoiiaioooaiiii 
+ aooifloioaioiaiio + aooiaoiiaiiooioo + aoioooiiaioiaioo 

+ +4(000000110101^110 + aooiaoiooiooaiii) (9) 

The polynomial invariants of four qubits were constructed in ( |Luque and Thibon 2003| ). 
Here, the Hilbert series is 

"(^) = (i-t^)(i-mi-t)a - (i«) 

For five qubits, the Hilbert series and a few fundamental invariants were obtained in 
ULuque and Thibon 2005| ). 



2.3. Covariants 

To construct the invariants, as well as for the more difficult problem of classifying the 
orbits, one needs of the classical notion of a covariant. A covariant $ of / is a multi- 
homogeneous GsLocc-invariant polynomial in the form coefficients a^^.i^ and in the 
original variables x- , that is, an invariant in some space 

$ e s'^'^\n) ® s"' {V*) «) ■ • ■ (g) s°"'{v*) , (11) 

where a is the multidegree of $ in the x^^ . 

Clearly, a covariant can be interpreted as an equivariant map «$ from the irreducible 
representation 

Sa{V):^S°''{V)'S)---(E)S°"'{V) (12) 

of G to S'^{T-i). Such a map is uniquely determined by the image of the highest weight 
vector Va of Sa{V). This highest weight vector is the coefficient of the highest monomial 
in $, classically called the source of the covariant. The coefficients of the other monomials 
form a basis of weight vectors in the image of m$ . 

The covariants form an algebra, which is naturally graded with respect to d and a. We 
denote by Cd-a the corresponding graded pieces. The knowledge of their dimensions Cd-a 
is equivalent to the decomposition of the character of S'^{TC) into irreducible characters of 
G, and the knowledge of a basis of Cd;a allows one to write down a Clebsch-Gordan series 
with respect to G for any polynomial in a. Also, it is known that the equations of any 
G-invariant closed subvariety of the projective space P(7i) are given by the simultaneous 
vanishing of the coefficients of some covariants. 

A modern introduction to classical invariant theory can be found in the book IjOlver 1999|l . 
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3. LUT-invariants from SLOCC-covariants 

3.1. General construction 

A generating set of the algebra of the polynomial covariants for the action of the SLOCC 
group can be in principle computed by a slight adaptation of the classical method (the 
Cayley Omega process, see, e.g., (jOlver 1999|l l. The covariants can be obtained recur- 
sively from the simplest one (the ground form /) 



/ = 



/ , '^il---ik^ii 



(1) 



Jk) 



(13) 



by iterating an operation called transvection, defined by 



(*,$)^i- 



= trr2'i,---r2 



(1) -^,(.)^(a;' 



/(I) 



,/W 



where 



x^{x- 



fix = det 



//(I) 



,x"('=)) 



)x 



dx'o 



a 

dx'-i 
d 



(14) 



(15) 



dx"o dx" \ 

and tr : a;', x" — > x. 

In practice, obtaining a description of the algebra in terms of generators and syzygies 
seems to be out of reach for more than four qubits IjBriand et al. 2003l|Luque and Thibon 2005| ). 
Nevertheless it may be always possible to compute the smallest covariants with relevant 
geometric properties. 

As already mentioned, a basis Covfe of the space of the polynomial SLOCC-covariants 
can be identified with a basis of highest weight vectors in the symmetric algebra S(7i), 
so that one can write 

0eCovfc 
where V^ denotes the irreducible representation of G whose highest weight vector corre- 
sponds to the covariant (p. 

Polynomial invariants under LUT (resp. LSUT) live in S{V'^'') S{V*'^'') and hence 
in 



(v.^v;,) 






(resp. {V^(»Vpy 

0,0'GCoVfe 



'), 



where degcj) denotes the degree of cf) in the variables ai-^,,,i^. 

Note that if ^ is a covariant whose multidegree in the auxiliary variables is (ni, 
the corresponding irreducible representation is 



V^^S" 



'-) (E) ■ ■ ■ (g) S''" (C^) 



If (p and (/)' are two polynomial covariants whose respective multidegrees are (rii 

and {mi, . . . , ruk), then V^ (8> V7/ contains polynomial invariants under LUT (resp. LSUT) 



(16) 
(17) 

(18) 

,"fc) 
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if and only if ni = mi ,■■■ ,nk — mk- Moreover, combining the previous abstract nonsense 
identifying covariants with G-highest weight vectors and G-equivariant maps, and the 
canonical antilinear isomorphism of a Hilbert space with its hermitian dual implies the 
following result: 

Proposition 3.1. Denoting by ^^i & basis of SLOCC covariants of degree d in the 
entries of the tensor and multidegree a in the auxiliary variables, one has: 

1 The scalar products {^dil^dj) "^ith respect to the auxihary variables, the ai^,,,ii, 
being regarded as scalars, form a basis of the space of LUT invariants. 

2 Similarly, the scalar products ($2il*^d' i) (where d' is not necessarily equal to d), form 
a basis of the space of LSUT invariants. 

The hermitian scalar product induced by the one of V can be calculated by the formula 

{xi ■■■Xm\yi ■■ -ym.) = perm((a;j|yj)) 
if {xi\yj) — 1 when Xi = yj and otherwise. 

This property should be of interest for the study of entanglement measures, which are 
special LSUT invariants. Indeed, expressing such a measure as a simple combination of 
scalar products of covariants with known geometric properties might lead to interesting 
insights. 

In the sequel, we will denote by CovsLocc(fc) (resp. Iuvlut, Ihvlsut) the algebra of 
polynomial SLOCC-covariants (resp. LUT-invariants, LSUT- invariants). Note that these 
algebras are multigraded. The space of multihomogeneous SLOCC-covariants (resp. LUT- 
invariants, LSUT -invariants) of degree n in the ai-^...i^, and d = (di, • ■ • , dk) in the aux- 
iliary variables (resp. degree n in the ai-^...i^'s and the ai^...i^'s, degree ni in the ai^.-.j^'s 
and degree 712 the ai-^...i^'s) will be denoted by CovsLocc(fc;"; d) (resp. InvLUT(fc;"), 

InVLSUT(^; ("-1,"2)))- 

3.2. Hilbert series 

From Proposition 13.11 we see that the knowledge of the Hilbert series of the SLOCC- 
convariants allows one to compute the Hilbert series of the LUT and LSUT-invariants. 
We will denote by 

hsLocc{k;z;u) = ^dimCovsLocc("'; ^; d)2:"u'^ (19) 

where u*^ — u^^ ■ ■ ■UjJ', the Hilbert series of InvsLocc- The Hilbert series of the algebras 
Iuvlut and Iuvlsut are obtained from /isLocc(fc; z] u) by the formulae 



n)z 



2n 



KuT{k\z) = 2^^dimInvLUT(fc;2 

= ft-sLocc(fc; z'^; u) hsLocc{k; z"^; u)|^_^ 

= C 1 z.tii....,Ufe l/lsLOCcifc; ^t; (Wl, . . . , Ufc))/lsLOCc('^i 7J (^1 5 ■ ■ • 5 ^fe )| J 

(20) 
where denotes the Hadamard product of the power series ring C [[z, ui, . . . , Uf.]] (that 
is, u" u'^ = i5q/3u"), and GTxi,...,x^f means the constant term of the series / with 
respect to the variables xi, . . . , a;„. 
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Similarly, one has 

KsvT{k]z) = X]„i,„2dimInvLUT(fc;("-i,?^2))^"'z"= 

= ^sLocc(fc;2;u)0c[[z,-]] /isLocc(fc;^;u)|^_^ 

(21) 
where 0c[[z,z]] denotes the Hadamard product in C [[z,z]][[ui, • • • ,Ufe]] (i.e. considering 
C [[z, z]] as the ring of scalars). 
Hence, 



dimInvLUT(fc;2n) = y^(dim Covslocc(?^; k] d))^ 



(22) 



and 



dimInvLsuT(fc; (n-i,"'2)) = y^ dim Covslqcc {ni]k\d) dim Covslocc ("-2 ; fc; d) . (23) 



Classical methods of invariant theory allows one to express the Hilbert series of algebras 
of covariants as a constant term (see IjBriand et al. 2003)l for an example). Hence, the 
Hilbert series of unitary and special unitary invariants are 



and 



n(i--r^) 



-2n2 



h^^rj.{k;z) = -jrCTt^u < 






11(1 - ury 



h-LsvTyk] z) — — ^Ci^ 



(24) 



(25) 



n [(i-zu")(i-zu") 

ae{-i,+i}'' 

These expressions have been first derived by Beth et al. (unpublished, see l|Grassl et al. 2002J) '). 
using a different method. 



4. Simplest invariants 

4.1. Dimension formulas for SLOCC-covariants 

The characters of the irreducible polynomial representations of the group GL(2, C )^ are 
the products 

sx := Sa(i) • • • SaC") (26) 

where A = (A*-^-*, . . . , A^*^^) is a tuple of partitions A*^*^ of length at most 2, and Sx(i} the 
corresponding irreducible character of GL{2, C ) i.e., a Schur function IjMacdonald 199l1l . 
In particular, the characters of the one-dimensional representations 

det'(g) = (det5(i))'Hdet5(2))'^ •■•(det5('=))''= , (27) 
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containing the SLOCC invariants, are the products 

Sillll)S(l2l2)---S{l^l^:) (28) 

and the character of GL{V) in S'^{V) is Sd- 

Hence, the dimension of the space of invariants of degree d and weight 1, which is also 
the multipUcity of the one-dimensional character det' in S'^(y), is given by the scalar 
product 

dimInVsLocc(rf,fc;l) = (srf | S(,j,j) ••■S(,,j^)) (29) 

of SLOCC characters (where 1 = (^i, . . . , Ik))- We see that InvsLocc(rf; k; 1) can be nonzero 
only if the condition 

d=2li = 2l2 = ■■■ = 2lk (30) 

is satisfied. Hence, 

dimInvs,occ(2/;A:)= ^ I^a'-'Xa, (31) 

Ah2! ^^ 
i(A)<2 

where x" denotes the value of the irreducible character x" (labelled by the partition 
{1,1)) of the symmetric group &21 on the conjugacy class A — (l™i2™2 • ■ • n™"), and 
zx = ni»'"'TO«!, cf. UMacdonald 199111 . 

In the same way, the SLOCC-covariants of a fc-qubit system form the algebra 

Cov = [S{V®^) ® S{V* © ■ • ■ ® T/*)]"^°^^^^ (32) 

which can be graded according both to the degree in the ai^...ik ^nd the multidegree in 
the auxiliary variables. A similar reasoning gives the dimension of the space of covariants 
of degree d as 

k 

dimCovs,occ(d;fc)=V-| y^x'l ■ (33) 



\i(A)<2 / 



Although impratical for finding closed forms of the Hilbert series, these expressions 
are useful for computing the first terms. 



4.2. Simplest SLOCC-covariants 

The space of covariants of degree 1 is generated by the ground form 

/= Y. a..-^.^---^- (34) 






The dimension of the space of covariants of degree 2 oi a k- qubit system follows from 
formula (|33|l 

dimCovsLocc(2,fc) = 2'=-i. (35) 

Observe that the only multihomogeneous covariants in this space have a multidegree in 
the auxiliary variables belonging to {0,2}'^. If d is any tuple, we will denote by |d|a 
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the number of occurrence of a in d. The dimension of the space of covariants of degree 
d = (di, • • • , dfc) S {0, 2}'^ in the auxihary variables is 

dimCovs.occ(2,fc;d) = ((x^)('="|d|o)(^n)|d|o|^2^ ^ | if jdjo is odd^ ^^^^ 

Hence, we can state the following result: 

Proposition 4.1. The space of the covariants of degree 2 of a fc-qubit system has di- 
mension 2'^'^^ and is spanned by /^ and the polynomials 

Bd = (/,/)^^--^^, (37) 

where d = (di, ■ • ■ , dk) G {0, 2}'^ and |d|o is even. 

Note that if k is odd, there are no invariants of degree 2 and if k is even the invariants 
are all proportional to the hyperdeterminant (/, fy^ ' . 

The dimension of the space of covariants of degree 3 is, from formula (|33|l . 

dimCovs^occ(3,A:) = ^3'=-^ + i. (38) 

The only covariants of degree 3 in the entries of the tensor have a multidegree d = 
{di, . . . ,dk) & {1,3}*^ in the auxihary variables. Let d = {di, . . . ,dk) G {1,3}'^ be a 
multidegree, the dimension of the space of the covariants having multidegree d is 

dimCovs,occ(3,fc;d) =((x'')l''l^(x')l''l^1(x')) .oQ^ 

= i(2ld|i-i + (_i)|dU) I'^yJ 

if |d|i > 0, and 

dimCovs,occ(3,fc; (3^=)) = ((x')'lx') = 1- (40) 

This implies that all the homogeneous covariants of multidegree (3") are proportional to 
/^. From H39|) . the dimension of the space of the /c-linear covariants of degree 3 is 

dimCovs,occ(3,fc;d) = i (2^-^ + (-l)'^') (41) 

Let us denote by {C'i}j^2,...,i(2'=-i+(-i)'=-M ^ basis of the space of covariants of multi- 
degree (l'"'). Applying transvections with the ground form to the Q, we obtain invariants 
of degree 4. Remark that the dimension of the space of invariants of degree 4 is equal 
to the dimension of the space of multilinear covariants of degree 3, so that we recover a 
result of ( |Brylinski and Brylinsky 2'002| ). We will denote by (Di) a basis of the space of 
SLOCC invariants of degree 4. 

4.3. Polynomial LUT-invariants of degree 4 

From ProDOsition l4.1l one can construct a basis of the space of LUT invariants of degree 
4. If d = (di, . . . , d/c), the dimension of Covslocc(2, fc; d) is or 1. Hence, the only 
possibilities are the squared norms 

Bd := (Sd|Bd) (42) 

where d = (di, • • • , dk) G {0, 2}^ and |d|o is even. 
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The dimension of the space is the coefficient of z'' in the Hilbert series 



dim InvLUT (4, k) 



Thus, the foUowing resuh holds: 



d 



(dimCovsLocc(2,fc;d))" = 2"- 



fe-i 



(43) 



Proposition 4.2. The space of the LUT invariants of degree 4 of a fc-qubit system has 
dimension 2''"^ and is spanned by the polynomials Bd and (/|/)^. 



Furthermore, one can show that 

B22...2 = (/'|/'>=2^-(/|/)2 



E 

d5^(2,...,2) 



B, 



(44) 



Indeed, (/|/)^ = if ^ f\f '^ /)j and the classical Clebsch-Gordan series allows one to 
express / (g) / in terms of transvectants, whence the result. 

If we consider only the LSUT group, one has generators of bidegree (4, 0), (3, 1), (2, 2), 
(1, 3), (0, 4) in the components of the state and their conjugates. The subspace of bidegree 
(2, 2) is the space of LUT-invariants of degree 4: 

InvLsuT((2, 2), /c) = InvLUT(4, k). (45) 

The subspace of bidegree (4, 0) if the space of the SLOCC invariants of degree 4 



Inv, 



,((4,0),fc)=InvsLocc(4,A:). 



(46) 



In the same way, the subspace of bidegree (0,4) is the space of the conjugates of the 
SLOCC invariants of degree 4. A basis of the space of the LSUT-invariants of bidegree 
(3, 1) can be obtained from the scalar products {Ci\f). 

Proposition 4.3. The subspace of LSUT-invariants of degree 4 of a /c-qubit system has 
dimension 

and is spanned by the polynomials Di (bidegree (4,0)), C^ = {Ci\f) (bidegree (3,1)), 
Bd (bidegree (2,2)), C^ (bidegree (1,3)) and A (bidegree (0,4)). 



5. Examples 

5.1. LUT-invariants and linear entropies 



Let 1^1') = J2 c-ii-.-ik 1*1 ■■ ■ ik) be a pure fc-qubit state. Meyer and Wallach ( |Meyer and Wallach 20021 ) 
have defined an entanglement measure Q by 



k 
i=l 



where 



D^\\^)) 



(ei, 



E 



fell*) 






(47) 
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In this expression, (ef |^) denotes the coefficient of |ei • • ■ ei-iSti+i ■ ■ ■ Ck) in |^), and the 
double bars mean the squared modulus of the determinant. 

The interest of this measure resides in its physical interpretation, related to the av- 
erage purity of the constituent qubits IjBrennen 2003|l or the linearised form of the Von 
Neumann entropy of a single qubit with the rest of the system. 

Emary has remarked ( |Emary 2004| ) that the functions Uj are entanglement mono- 
tones, thus in particular LU-invariants. Hence, each D]^ can be written in terms of 
squares of transvectants. We have 



7?« = 



2fc-2 



E 



Brf 



(48) 



d=(dl,-- ,dfc)e{0,2}'' 
d,=0 



Indeed, D]' = [|$i ® $i||^, where <&,; = i?22--202--2, and the result follows again from the 
Clebsch-Gordan series. 

Hence, in terms of our basis, the quantity Q(|^)) has the simple expression 



Qm) = 



1 



2^-'^k 



E 



IdLBd. 



(49) 



d=(di,---,dfc)e{o,2}'= 



5.2. LUT-invariants for 3-qubits 

The algebra of covariants of 3 qubits is generated by the polynomials fLe Paige 1881| ) 



Hx :— 



Hy-.^ 



H,:= 



t>yoOzo 
dyodzi 


dyidzo 
dyidzi 


dxndzQ 
dxodzi 


dxidzo 
dxidzi 


dxodyo 


dxidyo 



dxodyi dxidyi 



T := 



Of df 

dxQ dxi 
dxQ dxi 



A:=(T,/)^ 
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From these polynomials, we can construct the following LU- invariants 

Am := (/I/) 
B200 := {Hx\Hx) 

B020 '■— {Hy\Hy) 

B002 := {Hz\Hz) 

Cm :- (TIT) 

Dooo := (A|A) 

F222 := (A/'IT^) 

Grassl et al. l|Grassl et al. 2002}! have computed a minimal system of seven generators 
(denoted by fi) of the algebra of LU invariants. We shall give their expressions in terms 
of scalar products of covariants. 

The generator of degree 2, /i is clearly Am. To define generators of degree 4 and 6, 
the authors introduce the notation 

/cr,r,p := 2J aijkai-j-kp (50) 

i=(,l,>2,...,.„), 
J — {jl,J2''-'Jn), 



where 1" — 


(V(i),-- 


■ ,io{n)) and a^k = ai^j^ki ■ ■ •ai„j„fe„. Their gi 


and 6 are 




/2 :== /(12),(12),Id = Ajii - B200 - B020 
/s := /(i2),id,(i2) = Ajii - B200 - B002 
fi '■— /id, (12), (12) = Ajj^ — B020 — B002 



3 3 

h := /(I2),(23),(i3) = Am + 2^111 ^ 2^"i*-^200 

+B020 + B002) 

Note that these invariants appear in many places in the literature, for example in 
IKempe 1999| |. 

The generator of degree 8 is Dooo and the generator of degree 12 is 

/y:^ A([ll,00]{00,00}-[11,00]{11,11} 

+ [11, 01]{00, 01} + [11, 10]{00, 10} 

+2[11, 10]{01, 11} - 2[01, 00]{10, 00} 

-[01, 00]{11, 01} - [10, 00]{11, 10} 

-[10, 01]{00, 00} - [10, 01]{01, 01} 

+ [10,01]{10,10}+[10,Gl]{ll,ll}f 

where [1112, J1J2] = flj^j^oajiiai - aiii2iaiij20 and {11*2, J1J2} == ainaoajusi + «Jii2i«iii20- 
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Dn 



\GHZ) 


X 


X 


X 


X 


\W) 


X 


X 


X 





|Bi> = 1001) + 1010) 


X 











|B2> = |001> + |100> 





X 








IB3) = |010> + |100> 








X 





|000) 















Table 1. SLOCC orbits of three qubit states. 

With our notations, one has 

1 /3 

/? = -^Dooo ( X (B200 + B020 + B002) — Ajjj 

+2Cm — 4B200B020B002 + nF222- 

The authors of IjGrassl et al. 2002)1 obtained the Hilbert series using residue calcula- 
tions in Magma. We have been able to reproduce their results evaluating (|24|l using a 
very efficient algorithm due to Guoce Xin IjXin 2004|l in a Maple implementation. Sum- 
marizing, we have: 

Proposition 5.1. The algebra of local unitary invariants pure 3-qubit states is generated 
by Am, B200, B020, B002, Cm, Dqoo and F222- Its Hilbert series is 

l-t24 

^"""^^' '^ = (l-t2)(l_i4)3(l_i6)(l_i8)(l_il2) (51) 

which where the numerator reflects the existence of a unique syzygy in degree 24. 



5.3. Classification of the orbits under SLOCC transformations 

The normal forms of 3-qubit states under SLOCC transformations are known since 1881 
| |Le Paige 1881| |. As shown on tabled the SLOCC orbits can be characterized by the 
vanishing or non vanishing of a set of four LU-invariants. In the table, a x means the 
non-nullity of the invariant. Hence, (|48|l implies that in this case, the "onion classifica- 
tion" | |Miyake 2003| ) can be described only in terms of proper entanglement measures 
(entanglement monotones), see figH 



5.4. LS UT-invariants for 3-qubits 

Another unpublished result of Grassl et al. IjCrassl et al. 200'2|l can be recovered from 
H25|l by means of Xin's algorithm (jXin 2004|l . It is the Hilbert series of the algebra of 
LSUT-invariants of three qubits, 

z^z^ -I- z^z^ -I- z'^j'^ + 1 

(1 - ZZ)(1 - Z4)(l - z2z2)2(l - z4)(l - zi3)(l - z3z) ^ ' 
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Fig. 1. SLOCC orbit structure for 3-qubits 



«000 I 1001 I 1010 I 1011 I llOO I llOl I 1110 I iiii I 
3 + 3i I 3 + 3? I 3 + 3? | 2 + i | 3 + 2i | 1 + 2i | 2 + 3? | 3 + i | 



Table 2. Random values of the Qijk's 



This expression suggests that the algebra has a Cohen-Macaulay structure with 6 primary 
invariants of respective bidegrees ((1, 1), (0,4), (2, 2), (2, 2), (4, 0), (1, 3) and (3,1) and 3 
secondary invariants of bidegree (2, 2), (3, 3), (5, 5). The set of primary invariants is 



7' = {Aiii,/2,/3,A,A,s2 



, am, aooo with the random 



Computing the Jacobian of Am, /2, /a. A, A, S2, S2, aooo, 
numerical values given in table 12 one finds 

-53279560564736 - 243669580382208i ^ 0. 

This implies that the polynomials Am, /2, /s. A, A, S2 and S2 are algebraically indepen- 
dent. The set of secondary invariants is S = {f^, f^, /4/5}. The polynomial f^ (resp. /s) 
is linearly independent of all algebraic combinations of bidegree (2,2) (resp. (3,3)) of 
elements of V. Furthermore, one has two syzygies involving respectively /| and /|: 



8/1/5 - 6/4/2 + 3/1 - 3|A|2 + 3/1 - 6/4/3 + ft + 3/1 - 6/3/2 - 12|s2p = (53) 



and 



- 18/4/1^ - 18/3/f - 18/2/1^ + 11/1' + 18As^ - 36I52IV3 + 18As^2 _ 72/4/3/2 
+ 30/4/2/1' + 30/4/3/1' - 36|s2|'/2 + 60|s2p/i' + 3/I/1 + 3/I/1' + 30/3/2/1' - 36|s2|'/4 

+ 3/|/i'-3|Ap/i' + 16/1 = 0. (54) 

This implies the following property: 
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II ^ I 



II ^ I 



2 





4 





6 


6 


8 


46 


10 


110 


12 


344 


14 


844 


16 


2154 


18 


4606 


20 


9397 


22 


16848 


24 


28747 


26 


44580 


28 


65366 


30 


88036 


32 


111909 


34 


131368 


36 


145676 


38 


149860 


40 


145676 


42 


131368 


44 


111909 


46 


88036 


48 


65366 


50 


44580 


52 


28747 


54 


16848 


56 


9397 


58 


4606 


60 


2154 


62 


844 


64 


344 


66 


110 


68 


46 


70 


6 


72 





74 





76 


1 



















Table 3. Hilbert series of LUT invariants for 4 qubits: values of the a^. 

Proposition 5.2. The algebra of LSUT invariants of three qubits is a free module over 
a polynomial algebra (Cohen-Macaulay structure) 



Inv, 



^C[VV 



(55) 



ces 



5.5. LUT invariants of four qubits 

Again, we have computed the Hilbert series of LUT covariants of 4 qubits by means of 
Xin's algorithm. This allowed us to reproduce another result of l|Crassl et al. 20021 . 

P(z) 



r(4;z) 



Q{z) 



(56) 



with P{z) = 1 + ^- aiZ^z^ where the Oi are given in Table Inland 



Q{z) - (1 - z"')(i - ^Ta - n'a - ^ya - ^') • 

This suggests that the algebra has a Cohen-Macaulay structure with 19 primary invari- 
ants and 1449936 secondary invariants. The complete knowledge of the generators is with 
no doubt out of reach, nevertheless one can compute the first primary invariants using 
the covariants obtained in a previous paper IjBriand et al. 2003JI . The simplest one is the 
scalar square of the ground form 



Al 



111 



(/l/> 



There are 6 bi-quadratic linear covariants of degree 2 and 1 invariant. This allows to 
construct unitary invariants of degree 4: 



B2200 = 


(-B2200 -S2200 


B2020 = 


(-B2020 -S2020 


B2002 = 


(-B2002 -S2002 


B0220 = 


(-B0220 -B0220 


B0202 = 


(-B0202 -B0202 


B0022 = 


(-B0022 -B0022 


B = 


Soooo-Boooo 
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{i,i) 



(iJ) 



(iJ) I 



(iJ) I 



(iJ) 



(0,0) 


1 


(1,3) 


-1 


(2,2) 


2 


(2 


4) 


6 


(2 


6) 


9 


(2,8) 


4 


(2, 10) 


3 


(3,3) 


7 


(3,5) 


12 


(3,7) 


12 


(3,9) 


7 


(3,11) 


2 


(3, 13) 


-3 


(4,4) 


28 


(4,6) 


42 


(4,8) 


52 


(4, 10) 


36 


(4, 12) 


12 


(4, 16) 


1 


(5,5) 


43 


(5,7) 


79 


(5,9) 


92 


(5,11) 


36 


(5, 13) 


-1 


(5,15) 


-12 


(5, 17) 


-6 


(5,19) 


-1 


(6,6) 


132 


(6,8) 


199 


(6,10) 


161 


(6, 12) 


53 


(6, 14) 


-9 


(6, 16) 


-27 


(6,18) 


-10 


(7,7) 


214 


(7,9) 


236 


(7,11) 


129 


(7, 13) 


-12 


(7, 15) 


-83 


(7, 17) 


-63 


(7, 19) 


-15 


(7,21) 


-2 


(8,8) 


339 


(8,10) 


289 


(8,12) 


110 


(8, 14) 


-115 


(8,16) 


-169 


(8, 18) 


-82 


(8,20) 


-21 


(8,22) 


-3 


(9,9) 


306 


(9,11) 


160 


(9, 13) 


-154 


(9,15) 


-363 


(9, 17) 


-253 


(9, 19) 


-82 


(9,21) 


-12 


(9,23) 


3 


(10, 10) 


268 


(10,12) 


-96 


(10, 14) 


-513 


(10, 16) 


-510 


(10,18) 


-234 


(10 


20) 


-37 


(10 


22) 


12 


(10, 24) 


3 


(11,11) 


-126 


(11,13) 


-676 


(11 


15) 


-818 


(11 


17) 


-465 


(11,19) 


-85 


(11,21) 


76 


(11,23) 


41 


(11 


25) 


4 


(12 


12) 


-681 


(12, 14) 


-1045 


(12, 16) 


-763 


(12,18) 


-221 


(12 


20) 


133 


(12 


22) 


154 


(12, 24) 


36 


(12,26) 


3 


(13,13) 


-1152 


(13 


15) 


-985 


(13 


17) 


-359 


(13, 19) 


265 


(13,21) 


424 


(13,23) 


216 


(13 


25) 


41 


(13 


27) 


3 


(14, 14) 


-1094 


(14, 16) 


-543 


(14, 18) 


245 


(14 


20) 


705 


(14 


22) 


496 


(14, 24) 


154 


(14,26) 


12 


(14,28) 


-3 


(15 


15) 


-569 


(15 


17) 


318 


(15, 19) 


1058 


(15,21) 


992 


(15,23) 


424 


(15 


25) 


76 


(15 


27) 


-12 


(15, 29) 


-2 


(16, 16) 


233 


(16,18) 


1188 


(16 


20) 


1334 


(16 


22) 


705 


(16, 24) 


133 


(16, 26) 


-37 


(16,28) 


-21 


(17 


17) 


1333 


(17 


19) 


1734 


(17,21) 


1058 


(17,23) 


265 


(17,25) 


-85 


(17 


27) 


-82 


(17 


29) 


-15 


(17,31) 


-1 


(18, 18) 


1736 


(18,20) 


1188 


(18 


22) 


245 


(18 


24) 


-221 


(18, 26) 


-234 


(18,28) 


-82 


(18,30) 


-10 


(19 


19) 


1333 


(19 


21) 


318 


(19, 23) 


-359 


(19,25) 


-465 


(19,27) 


-253 


(19 


29) 


-63 


(19 


31) 


-6 


(20, 20) 


233 


(20, 22) 


-543 


(20, 24) 


-763 


(20 


26) 


-510 


(20 


28) 


-169 


(20, 30) 


-27 


(20, 32) 


1 


(21,21) 


-569 


(21 


23) 


-985 


(21 


25) 


-818 


(21, 27) 


-363 


(21,29) 


-83 


(21,31) 


-12 


(22 


22) 


-1094 


(22 


24) 


-1045 


(22, 26) 


-513 


(22, 28) 


-115 


(22,30) 


-9 


(23 


23) 


-1152 


(23 


25) 


-676 


(23, 27) 


-154 


(23, 29) 


-12 


(23,31) 


-1 


(23 


33) 


-3 


(24 


24) 


-681 


(24, 26) 


-96 


(24, 28) 


110 


(24, 30) 


53 


(24 


32) 


12 


(25 


25) 


-126 


(25, 27) 


160 


(25, 29) 


129 


(25,31) 


36 


(25 


33) 


2 


(26 


26) 


268 


(26, 28) 


289 


(26, 30) 


161 


(26,32) 


36 


(26 


34) 


3 


(27 


27) 


306 


(27, 29) 


236 


(27,31) 


92 


(27,33) 


7 


(28 


28) 


339 


(28 


30) 


199 


(28, 32) 


52 


(28, 34) 


4 


(29, 29) 


214 


(29 


31) 


79 


(29 


33) 


12 


(30, 30) 


132 


(30, 32) 


42 


(30, 34) 


9 


(31 


31) 


43 


(31 


33) 


12 


(32, 32) 


28 


(32, 34) 


6 


(33,33) 


7 


(33 


35) 


-1 


(34 


34) 


2 


(36, 36) 


1 























Table 4. Hilbert series of LSU invariants for 4-qubits: values of the a^ 



The polynomial (/^|/^) is algebraically dependent of the other ones: 

(/ 1/ ) — 16A — (B2200 + B2020 + B2002 + B0220 + B0220 + B0202 + B0022 + B) . 
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The space of linear covariaiits of degree 3 is spanned by two quadrilinear polynomials 

L-iiii - U1-D2020J 
and four cubico-trilinear covariants IjBriand et al. 2003JI 

n — { f n ^0100 

(-^3111 — U,-O2200j 

1000 



C13II — if, -B2200) 
*-^1131 — Uj -02020^ 

r — If n ^looo 

L'1113 — \J, i'2002j 

With these polynomials, one can construct a set of twenty generators for the space of 
unitary invariants of degree 6: 

A3, 

AB, AB2200; AB2020: AB2002, ABo2205 AB0202, AB0022, 
(CiiiijC'iiii), (CiiiijC'iiii): (Ciiii,/-Soooo), (C'iiii,Ciiii), 
(^iiiu ^1111)1 (C11111/-S0000), (/-Soooo, C'liii), {/^oooojC'iiii}, 
(C'31111 C'3111), (C1311, C1311), (C1131, Ciiai), (C1113, C1113). 

The series suggests that the algebra has a Cohen-Macaulay structure with 19 primary 
invariants (one of degree 2, seven of degree 4, four of degree 8 and one of degree 10). The 
polynomials 

Aim, 

B, B22001 B20205 B2002, B02201 B0202, B0022 

(CiiiijC'im}, {Ciiii,AB), (C3iii,C3iii), (Ci3ii,Ci3ii), (Cii3i,Cii3i}, (Ciii3,Ciii3} 
(-D4000, -D4000), (-D040O) ^040o), (-D004O) -D0040), (-D0004, -D0OO4) 

(-E'311l7-£'311l) 

where 

^4000 — (^, C'3111) , 

-^0400 = (^7^1311) , 

£'004Q ~ {A, C1131) , 

£'0004 = {A, C1113) , 
-D2200 = (^J C'3111) 

?220o) 

are algebraically independent and hence good candidates to be primary invariants. 



and^3iii-(A^2---^""° 



5.6. LSUT invariants of 4 qubits 

Finally, we can compute the Hilbert series of LSUT invariants of 4 qubits by the same 
method, and again recover a result of IjGrassl et al. 2002}! . 

/ilsut(A:;z,z) = — — (57) 
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with P{t) ~ ^. aijZ^'z-' , the a^ being given in Tableland 

Q(i) - (1 - Z-Z)(\ - 02^2)4(1 - zH'^^l - Z^){\ - z4)2(l _ z6) 
(1 - z2)(l - Z*)2(l - z6)(l - z3z)3(l - Z-Z^f{\ - z"^) 
{\ - Z'^T^){\ - Z-Z^){:\ - Z^^)) 

6. Conclusion 

We have proposed a new method to compute bases of the algebras of unitary invariants 
of qubit systems. This method involves as an intermediate step the calculation of the 
SLOCC covariants, which have a more transparent geometrical meaning (at least in small 
degrees), and leads naturally to new bases in which the known entanglement measures 
tend to admit rather simple expressions. 

The complete description of the algebra of unitary invariants for pure /c-qubits is 
definitely out of reach of any computer system for fc > 3. This impossibility means that 
such a study is not physically relevant and that only a few invariants with interesting 
geometrical properties will be significant in the realm of quantum information theory. 
Finally, a natural question is whether these constructions can be extended to mixed 
states. 
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